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The dynamics of systems undergoing a second order 
symmetry breaking phase transition has been studied re- 
cently in several instances particularly in associa- 
tion with the theory of topological defect formation p|| 
in cosmology (t]] and in experiments in 3 He || and ^ie 
H . With experimental data in the pipeline for the colli- 
sions of heavy ions at RHIC and later at CERN the effort 
to identify phase transitions in nuclear matter is also a 
very active field of research. The Langevin dynamics of 
effective scalar theories (such as the a model) affords one 
of the few quantitative windows into such transitions p0[ . 

Theoretically this situation has been modeled by the 
out-of equilibrium dynamics of a classical scalar field the- 
ories (with a given number of flavors) in contact with an 
external environment at a temperature T, which may be 
a function of time. The environment is only known sta- 
tistically and its behavior can then be described in terms 
of stochastic fields obeying a fluctuation-dissipation re- 
lation. The effective evolution of the fields is therefore 
described by a Langevin field equation. 

This stochastic classical description is of course only 
an approximation to the full quantum evolution. Re- 
cently attempts to treat quantum field theories dynam- 
ically jll|] have been made but the possibilities afforded 
by these approximations thus far do not allow for the 
correct description of thermalization processes |12 (with 
the exception of Boltzmann transport theory |13|) and 
have therefore very limited applicability to the descrip- 
tion of the long time dynamics of a system undergoing a 
symmetry breaking phase transition. 

In addition a reasonable case can be made that close 
to a second order transition the theory is effectively at 
high temperature T >> m(T), where m(T) is the tem- 
perature dependent mass scale, and its infrared should 
behave approximately classically. In this regime some 
quantum effects can also be included effectively through 
the appropriate choice of the coupling of the fields to the 
external environment (the quantum fluctuations in this 



case) and/or their contributions to the mass and cou- 
plings of the dynamical fields. This line of considerations 
|14f has been used recently to describe the effective evo- 
lution of the long wave-length modes of the color fields 
of non-Abelian gauge theories in the high temperature 
regime, as a stochastic classical field theory. 

In cosmology and high energy physics field theories 
are relativistic. This implies that the time evolution of 
scalar fields is governed by a second order derivative in 
time. The canonical Langevin equations however are usu- 
ally formulated in terms of a first derivative evolution in 
time. This difference is unimportant in equilibrium, but 
modifies real-time properties. In certain circumstances, 
which I discuss in detail below, the second order evo- 
lution effectively leads to a first order (or overdamped) 
description of certain parts of the system. 

The general properties of Langevin field evolutions are 
that, under very minimal conditions, the fields equili- 
brate for long times to the canonical (classical) Boltz- 
mann distribution. The dynamics are very rich, realizing 
and generalizing all the near canonical equilibrium prop- 
erties of the classical field evolution, including in partic- 
ular the features captured by the perhaps more familiar 
time-dependent Ginzburg Landau formulation (TDGL). 

In this paper I present a didactic introduction to 
the underlying formulation of a second time derivative 
Langevin system and its associated Fokker-Planck equa- 
tion. Its general dynamical properties are compared to 
those of other important classes of dynamics such as 
TDGL and the stochastic non-linear Schrodinger equa- 
tion (sNLS). This is done mostly in section |. (The treat- 
ment of the sNLS system is given in Appendix A.) In sec- 
tion [n] I derive analytic solutions of the Fokker-Planck 
equation in the particular case of the harmonic poten- 
tial. These solutions are fairly simple but to the best of 
my knowledge have not been discussed in the literature 
of stochastic field theories. A mean-field formulation of 



the dynamical problem is given in section III and it is 
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found that its solution is similar in effort to solving the 
associated Langevin mean-field equation, but with the 
stochastic averages already taken into account. In sec- 
tion |^ I discuss applications of the solutions of section 
[H| to the understanding of the full non-linear dynamics of 
the system. Although qualitative it will be argued that 
these solutions capture the essentials of the full dynam- 
ics in the critical domain and permit an understanding of 
the relevant thermalization times. A brief discussion of 
spinodal instabilities is also included. Finally in section 
^ I summarize the conclusions of this paper. 



I. THE LANGEVIN AND FOKKER-PLANCK 
EQUATIONS FOR SCALAR FIELDS 

We start by formulating the problem. The second or- 
der Langevin equation for a scalar fields </>(x, t) with an 
arbitrary interaction potential V((j)) is given by 



(d 2 t - V 2 ) d>(x) + + #(x) = £(x, t). 



(1) 



The stochastic fields £(x, t) are taken to be Gaussian and 
white, characterized by 

(£(x, t)) = 0, (£(x, f)£(x', 0) = fW 3 (x - x>(* - t'). ( 2 ) 

Although we have written this equation for a single scalar 
field <j>(x, t) the generalization of Eq. and of what 

follows, to a 0{N) symmetric theory, with N the number 
of flavors, is straightforward. 

We can proceed to reduce the order of this differential 
system, by introducing field generalized momenta conju- 
gate to tj>: 



d t 4>(x) = n(x), 

d t Tr(x) = -ipr(x) + V 2 c/>(x) 



SV{<t>) 
84>(x) 



(3) 
(4) 



An interesting limit of Eqs. arises when rj is large, 
or more precisely when dtit{x) << t]tt(x). Then we can 
write 



V d t( f>(x) = V 2 <t>(x) 



5<j){x) 



(5) 



Although this is the more conventional form of the 
Langevin field equation |l50] I will refer to it as the 
overdamped limit of Eqs. (£j]4|). 

The objective of the Langevin analysis is to measure 
the expectation value (over the stochastic fields) of any 
functional p of the fields n(x), <j>(x), generated by the evo- 
lution of Eq dlH). With the choice Eq. (|) this average 
is a Gaussian functional integral of the form 



(6) 



where the £ subscripts in the field and its momentum de- 
note their functional dependence on the stochastic fields 
t), through the evolution of Eq. (|!-[|). Relation (||) 
and (|!) can be generalized to fields £ with more compli- 
cated Gaussian distributions. 

In the presence of the stochastic fields £(x, t) the initial 
conditions for the fields 7r(x), 4>{x) may only be know 
statistically themselves. They can then be expressed in 
terms of a functional probability distribution P[tt, tfi], at 
the initial time. 

Throughout the evolution we assume that a time 
dependent probability distribution Ppp[(j)(x),Tr(x),t], a 
functional of the time independent fields and a function 
of time, exists. This is the Fokker-Planck probability 
distribution, which we require, as usual, to be positive 
definite and normalizable in the sense: 



Af / DirDcj)P F p{(t),iT,t} = 1. 



(7) 



Note that in general Af as well as other expectation values 
of the fields may be formally infinite as the ultraviolet 
cutoff of the theory is taken to zero. 

The time dependent expectation value over the 
stochastic fields £ of any quantity p[<fi, tt] is then given 
by 



(p){t)=M / DnD(/>P F p[(l>,n,t]p[<l>,ir] 



(8) 



In order to be useful this picture requires the explicit 
knowledge of Ppp[(j),7r,t]. The Langevin equation gen- 
erates a time-dependent (Fokker-Planck) probability dis- 
tribution Ppp[n, (j), t], which we can write formally 



Ppp[ir,<f>,t] = 

d D x8[fr(x,t) — tt(x)]S 4>(x,t) 



x) ), (9) 



where the brackets denote, as before, average over the 
stochastic fields. The static fields are the arguments of 
P F p[ir, 4>, t], whereas n(x,t), (f>(x,t) obey the Langevin 
Eqs. (§-(§. 

The probability P F p[ir, (j)^], obeys a functional 
Fokker-Planck evolution equation that can be computed 
directly by differentiating Eq. ([)]). Following Zinn- Justin 



dtPpp = 



s s 

d D x d t (f)(x,t)— h d t Tr{x, t)- 



5(j)(x,t) 1 " 1 ' '6%(x,t) 
X6[n(x,t) -n(x)]S ${x , t) - (j)(x)}) . (10) 



The properties of the 8 allow us to trade the functional 
derivatives of Tr(x,t),^)(x,t) for derivatives relative to 
tt(x), <j>(x). These can be taken out of the stochastic av- 
erage (...). The only remaining term inside this average 
results from the appearance of the stochastic fields £ in 
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the equation of motion for Tt(x,t). This term is of the 
form 



d u x£(x,t)5 [ir(x,t) - tt(x)]). 



(11) 



The Gaussianity of the stochastic fields enforces the iden- 
tity H 



(£(x,t)6[Tr(x,t)-n(x)}) = Cl( 



6[w(x,t) - n(x)]) 



(12) 



' <5£(x, t) 7r(x, t) 

Once again we can trade the functional derivative relative 
to Tr(x,t) for another relative to n(x). The expectation 
value of w^rr can in turn be obtained from the formal 
integration of the equation of motion by using a regular- 
ized delta function in time. We must effectively take half 
of this delta function to obtain 

5%(x,t) _ 1 
2' 



(13) 



Bringing together all the terms results finally in 

d t P FP [iT,^t] = -H FP P F p[ir,<f>,t], (14) 

where 



6n 



7/7T - V 2 



T) < 15 » 



It is useful to think in terms of the functional operator 
Tipp as the generator of infinitesimal time displacements 
of the probability functional Ppp . If, as in most applica- 
tions, the potential V((f>) is explicitly time independent 
we can invoke a separation ansatz for Pfp such that 



P FP [ir,<j>,t]=V[n,<j>}T(t) 



(16) 



Thus we can regard Eq. (|1J) as an analog of a functional 
Schrodinger equation, in imaginary time. Then we can 
write the time independent and dependent equations 



H F pP n = E n V n , d t T(t) = -E n T(t). 



(17) 



The functional dependence on the fields is now limited 
to the static probability eigenfunctionals V n . The time 
evolution of the Fokker-Planck distribution is completely 
characterized by the spectrum of eigenvalues of 7Yfp, E n . 
An orthogonal complete basis of functionals B n can in 
general be constructed from the set P n Q Note that it is 



*The set of functionals P n will only be guaranteed to be an 
orthogonal basis if Hfp is a Hermitian functional operator, 
which is not true in general. The relation between the B n and 
P n is usually very simple involving a power of the canonical 
distribution and it follows that E n are also the eigenvalues of 
the set B n . See eg. Jl^ for some examples in few degree of 
freedom problems. 



not possible in general to separate Eq. ( |l6| ) further into 
two equations, one in it and another in cj> via an ansatz 
like P F p = R[n]S[4>]. 

On general grounds we expect a steady state corre- 
sponding to thermal equilibrium to be reached for long 
times. Formally, we can then project the evolution of 
Pfp in terms of the eigenvalues E n and functionals B n 
as: 



P FP [ir,<i>,t] = ^C n B n [ 7 



-E„t 



(18) 



n=0 



where the C;'s are the projections of Pfp at the initial 
time onto the basis of eigenfunctionals B n : 

C n = J DvD(f)B n [TT, 4>]P FP {Tr, <f>,t = 0], (19) 

where both B n and Pfp[tt, 4> 1 1] are taken to be normal- 
ized. 

If equilibrium is approached for long times the corre- 
sponding distribution must be associated with the zero 
mode of Hfp, E — i.e. it is stable. Contributions from 
eigenstates V n , re^O, vanish exponentially, in a charac- 
teristic time i e q ~ E~ x , provided that the real part of 
E n is positive. This follows from the requirement that 
V(4>) is bounded from below. For the E n , solutions of 
the second time derivative eigenproblem there may be in 
addition an imaginary part. Thus the excited states de- 
cay away in time as an exponentially damped oscillator. 

The eigenfunctional Vq, characterizing the long time 
field probability is then interpreted as the equilibrium 
distribution. Taking Eq — in Eq. (17) we find 



■ A/"exp 



-0 



ir(x) 



(20) 



with S(<f>) = (i(V0(a;)) 2 +V(4>)). Here we took CI = 
2rj//3, which is the analog of Einstein's relation for Brow- 
nian motion, and ensures the balance between fluctu- 
ations (sourced by the fields £(x, t)) and dissipation 
(through the —i]ir(x,t) in Eq. (Q)) at long times. 

It is interesting to note that it can be read directly 
from Eq. (EOh that the momentum equilibrium variance 



(tt(x) 2 ) = 1/0 = T, 



(21) 



which expresses equipartition in a relativistic classical 
field theory. This quantity is also an excellent thermome- 
ter for the dynamical evolution. 

Now, because the canonical distribution is Gaussian 
in 7T we can proceed to obtain a reduced distribution in 
terms of (f> alone. Performing the Gaussian integral over 
all 7T we obtain 



P eq [<f>] =M'exp 



-0 I d D x S(<j>) 



(22) 



which is the canonical Boltzmann distribution. P e g[0] 
could have been equally obtained from the static solu- 
tion of the Fokker-Planck equation associated with the 
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perhaps more usual first derivative Langevin equation 
(see below). In this respect we see that in equilibrium 
it makes no difference to perform the evolution, with or 
without the second time derivative in Eq. (Q). Away from 
equilibrium, or even if we simply wish to study the fluctu- 
ations around it, we should take the appropriate form of 
the evolution, choosing to keep or neglect the second time 
derivative according to the physical picture under consid- 
eration. The difference between both of these evolutions, 
is expressed in terms of the higher eigenvalues and eigen- 
functionals associated with the different Hamiltonians, 
which I discuss below. 

For a harmonic potential V[<j>] — ^(fi 2 we can also 
read the equilibrium thermal propagator directly from 
Eq. (p2f). It is, in momentum space 



T 



P 



171 



(23) 



This is the canonical free Boltzmann propagator, which 
is the correct description for ideal classical waves at finite 
temperature. 

It is illuminating to compare the above framework di- 
rectly with that for the overdamped Langevin evolution. 
Following the same, but somewhat simpler, procedure it 
can be shown |16| that the Fokkcr-Planck Hamiltonian 
takes the form 



7~Lfp = — tt 



—A. v 2 + — 

2r\ 5(f> 8<f) 



loading to the eigenvalue functional equation 

A- — 

rid t P[4>]=-E N P[4>}. 



n 6 „ 2l 5V 



E N \ P[4>] = 0, 



(24) 

(25) 
(26) 



Again the canonical distribution is the eigenfunctional Pq 
corresponding to a zero eigenvalue, i.e. 

E = 0; P eq [<j)]=M'e-^ dD ^ v ^ +v ^\ (27) 

where as before the Einstein relation Q — 2r)/[3 must 
hold. 



II. APPROACH TO EQUILIBRIUM 

A complete dynamical solution of the ensemble (un- 
der average over the noise) of fields obeying the Langevin 
equation (^-^) can be obtained in principle by solving the 
corresponding Fokker-Planck equation. Unfortunately 
the task of solving for these non-linear functional equa- 
tions in general is quite monumental. The main difficulty 
is connected with the non-local form of the equations: in 
x-space this arises from the Laplacian term, while in k 
space the difficulty is transfered to the non-linear terms. 

Below I give closed form solutions to the equations in 
the fields and their conjugate momenta for the particular 



case of the harmonic potential. These are still interesting 
since they offer an alternative to the direct solutions of 
the Langevin equation, which, clearly, can also be found 
in the linear case. In the former picture, however, the 
averages over the stochastic fields are already taken into 
account. 

It is simpler to build some intuition for the solution of 
the overdamped Fokker-Planck equation first to which I 
now turn. 



A. The overdamped Fokker-Plank equation and its 
solutions 

The advantage of the harmonic case is that in Fourier 
space different modes decouple in the Langevin equation: 



where the stochastic fields obey 

(&(*)> = o, Mt)Z-k(t')) = ns(t-t') 

and 

d D k 



{x,t) = 



(27T) 



D 



(28) 
(29) 

(30) 
and simi- 



The fact that <p(x, t) is real implies 4> k = <fi*_ k 
larly for Because of these relations it is more conve 
nient to work with <f>^ = Ke((f> k ) and (\>\. — lm(0fe), which 
obey the Langevin Eq. (Eq) with £ 



R,I 



(&(*) B>I > = o, (C \m' I (t')) = "js(t-t'). (31) 

The Fokker-Planck equation for the modes now is 

r)d t P k = -H k P k (32) 



,R,I 



4r?, 



+ (k 2 + m 2 ) <f> 



R.I 



(33) 



which holds, as indicated, both for the real and imaginary 
components of <j> k . The solution of the Fokkcr-Planck 
equation then follows as 



The time independent solution is clearly 



P fe (^)P fc (^)ocexp 
which leads to 



-f3(k 2 



(34) 
(35) 



Pa = 



Afexp 
Afexp 
A/"exp 



-P 



d D k 

50 d D k 



(k + m 
k 2 ^ 



(36) 



D ^ 

1 ,„ , / „2 . II) 



{2k) d r 2 
d D x £ {V<t>{x)Y + 
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To compute the eigenfunctions corresponding to non- 
zero eigenvalues we need to solve the eigenvalue problem 



1 



2/3, 



,R,I* 



(k' + m 2 )^' 1 - 



— + En 



F 



N 



(37) 



where we wrote = PqFn- This equation has a familiar 
solution. To see this explicitly we perform a change of 
variables to obtain 



SX 



2X- 



2E 



N 



5X 



N 



(38) 



where X = y/ ' (3{k 2 + m 2 )^ 1 . Eq. (§|) has the solution 



E N = N(k 2 + to 2 ) 

F N = H N (X) = H N (y/0(tf + rf)4t*> 1 ) 



(39) 
(40) 



where Hn is the Hermite polynomial of order N . 

From the properties of the Hermite polynomials we 
know that these solutions are orthogonal under the mea- 
sure e~ x , which is just the k-part of the canonical distri- 
bution. The set of functions {Hn(X)c~ x Z 2 } thus con- 
stitutes an orthogonal basis. It can additionally be shown 
trivially that this basis is complete. 

The non-zero Eigenvalues Em set the time scales for 
the approach to thermal equilibrium, which we will dis- 
cuss in more detail below. 



B. The second order Fokker-Planck equation 

For the second time-derivative harmonic evolution the 
Langevin equation becomes 



We can bring Eq. ( p4| ) to the form of Eq. (|3 
the change of variables X± = a±7rv ' + b±<, 



|) by making 



with 



2 
a ± = O 



l±Wl-4 



k 2 + to 2 



7,2 2/3 ,. 2 . 2\2 / 

b ± = —{k + to ) I 



1±4 1-4 



k 2 + m 2 
j] 2 



The solutions of Eq. (44) then are 



E N 



P 2 



1± 



k 2 



(45) 



(46) 



(47) 



F 



N 



H N [X-i 



These solutions can be brought to products of Hermite 
polynomials in ir k ' and cf> k ' alone through the use of the 
property: 

N / i \ 

H N [x + y] = 2~ N / 2 J2 [NJ H N -i[x\/2}Hi[yV2]. (48) 



This property allow us to see that the solutions 



Pn = 



(F N [X+]t 



-E+t 



F N [XJ\e- E »*)P 



(49) 



form a complete orthogonal basis. 

A few properties of these solutions are worth point- 
ing out. The coefficients a±,b± and the eigenvalues 
E N can now be complex numbers. This happens for 
4(fc 2 + to 2 ) > ?f. In this case X+ = X*_ , E+ = {E N )* 
and the combination in ( |49"| ) remains real as required. 

The overdamped limit is recovered as j; 2 > > k 2 + to 2 . 
Then 



Cl(pk = TTk 
dlT k = -?/7Tfc 



(A: 2 + m 2 )0 fc +^ fe , 



(41) 
(42) 



which as in subsection II A can be written in terms of real 



and imaginary components. The Fokker-Planck Hamil- 
tonian for the latter is 



Hk — — 



n s 2 



R,I 



5 



4 s*F J 



(43) 



Sty 



R,I 



To find the excited states we proceed as before to factor 
out the canonical distribution P = F^Pq to get 




E-=N- 



V 



(50) 



(3 k 2 + to 2 
2 rf 



1 " «6 2 _ = f(fc 2 +TO 2 ). (51) 



The eigenfunctionals associated with ~ Nr][l + (k 2 + 
m 2 )/r) 2 ] are, in contrast, damped away rapidly. 

The characteristic thermalization times of the system 
are thus set by£jv = V-^jv- For the long wave length 
modes (k ~ 0) we have that the equilibration time t eq ~ 
r\j(k 2 +to 2 ). In the converse limit, of short wavelengths, 
the decay of all states is controlled at leading order by 
t eq ~ 2/to which is scale independent. In this sense the 
evolution of the short wavelength modes is quantitatively 
different from the overdamped Langevin system. 



III. A MEAN-FIELD APPROXIMATION TO THE 
NON-LINEAR FOKKER-PLANCK DYNAMICS 

It would no doubt be desirable to be able to solve the 
Fokker-Planck equations for more general (interacting) 
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potentials. At the next level of complexity it is possi- 
ble to attempt a solution of the Fokker-Planck equation 
in a mean-field approximation, which we discuss in this 
section. 

In order to illustrate the procedure let us consider 
an interaction potential of the form Vm — j(j)(x) 4 . In 
momentum space this potential generates a term in the 
Fokker-Planck equation for (f>k of the form 



= A 



d D 



P 



d D c 



{2tt) d (2 7 r) i 



p p q> q q>k- P -q- 



(52) 



A mean-field approximation to this term can be obtained 
by taking p = —q and the noise average such that 



SVmf 



= 3A 



d D 



P 



fa = 3XG(t)4> k . (53) 



Under this approximation the evolution equation can be 
written in terms of a mean-field Fokker-Planck operator 
where Vmf substitutes the potential V, with 



Q \ 

V MF = Vo + — G(t) 2 , 



(54) 



where Vo corresponds to the purely harmonic case, to- 
gether with the self-consistency condition 



G(t) 



P 



(27T) 



D 



d(j)p (j)p4'-pP[' 7T q, 4>p, t] 



(55) 



It is then clear that the under this Gaussian approxi- 
mation the theory remains harmonic, but with a time- 
dependent (self-consistently determined) mass. 

Unfortunately the time dependence of the mean field 
potential destroys the separability of the solution into 
a function of time, and a functional of the fields and 
their conjugate momenta. Closed form solutions are thus 
difficult to construct, which is of course also the case for 
the Langevin equation under a similar approximation. 
Numerical solutions of the this mean field Fokker Planck 
equation can nevertheless be easily obtained. 

As a starting point let us again consider the over- 
damped case first. Mean field approximations correspond 
in general to Gaussian probability distributions. In the 
overdamped case, this can only lead to 



Pmf [4>k ,t]=Af k (t) exp 



(56) 



where the normalization Af(t) is a function of A 2 (t) 

JV fc (t) = A fc (t)/VSF, 



(57) 



and must therefore be time-dependent. This does not of 
course introduce any additional dynamical freedom. 

The mean-field Fokker-Planck equation for PMF[4>k,t] 
translates into an equation for the coefficients Ak(t) of 
the form 



dAk r 2 
G(t)-- 



3XG(t) - A 2 T] A k , 



d D p 1 
j2^Aj- 



(58) 
(59) 



Clearly solving the Fokker-Planck equation in this mean- 
field approximation requires the solution of a system 
of coupled ordinary differential equations, one for every 
mode k. This is similar in effort to the solution of the 
Langevin mean-field equation with the advantage that 
the stochastic average is already taken into account. This 
latter step corresponds to the functional integral neces- 
sary to compute G(t) in the Fokker-Planck picture, which 
can be given in closed form because Pmf is Gaussian. 

The canonical (mean-field) distribution is of course a 
static solution of Eqs (p8-b3) where Ak obeys 



4 = 



T 



4 



-lD 



P 



1 



(60) 



It is then clear that for vanishing interactions A — 0, the 
free Boltzmann propagator is the solution of Eq. (^) . In 
general because the mean-field contribution is momen- 
tum independent Eq. ( p0| ) are equivalent to the 'gap' 
equation 



Am 2 = 3A 



Am 2 



T 

d D p 



T 



(2ir) D k 2 +m 2 + Am 2 ' 



(61) 
(62) 



For the second order Fokker-Planck equation the most 
general Gaussian function of the field and momentum 
modes is of the form 



PMF[^k,(t>k,t] = 

I A 2 

A4(i)exp ' 



(63) 



B 2 

Y<t>k4>-k + -^7T fe 7r_ fc + ClRe(TT k (t)-k) 



The normalization is now a function of all A 2 , B 2 , C 2 



N k = -JA 2 ,B 2 



^k-^k 



(64) 



The equations of motion for A\,B\ and C\ then follow 
from the second order Fokker-Planck equation: 



dAk 
dt 

dB k 
dt 

dC k 
dt 



where 



Ck 
A k 



k 



2 i 2 

1 m 



3AG(i) - rjTCl] 



= V B k [1 - TBI] 



[1 - 2TB 2 ] 



2Ck 



G(t) 



d D p 



3XG(t)) B 



I] . 



(2n) D A 2 p B 2 - C± p 



(65) 
(66) 

(67) 

(68) 
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This system has as static solutions the mean-field canon- 
ical distribution namely 



Al = 



k 2 



Am 2 (T) 



T 



Bl 



Cl 



l/T 
0, 



(69) 

(70) 
(71) 



where the equation for Ak is simply (|60[) . Other static 
solutions exist though. Additional stationary distribu- 
tions distinct from thermal equilibrium also occur in the 
microcanonical evolution of the class of models consid- 
ered here fl^| and seemingly result from the truncation 
of the system at any finite order in a hierarchy of corre- 
lators. The overdamped limit, k 2 + m 2 + Am 2 (i) << r] 
is obtained when C 2 ~ A\jr\ «B\,A\. 

Again we see that the effort in solving the mean-field 
second order Fokker-Plank equation is comparable to its 
Langevin counterpart, but with the stochastic averages 
taken into account. A renormalization scheme for ren- 
dering these equations ultraviolet cutoff independent can 
be constructed, if desired, following standard techniques 
for mean- field evolutions [fllT| . 



A is the ultraviolet cutoff which in many cases has phys- 
ical meaning because generally scalar field theories are 
effective low energy models. An example is the scale of 
separation between the true quantum behavior of excita- 
tions at large k 2 and the effective classical dynamics for 
the long wave-length modes which we describe stochas- 
tically. In this case, at high temperature and small cou- 
pling A, A ~ T. 

The interesting fact about m 2 (T) — > at criticality 
is that it supplies us with an arbitrarily long time scale 
separation between the thcrmalization time for the short 
and long wave-length modes. As we discussed in section 
H the long wave length modes k 2 ~ thermalize in a 
typical time-scale 



V 



m 2 (T) 



■ 00, 



(73) 



which is the expression of critical slowing down in our 
stochastic system. On the other hand the short-wave 
length modes thermalize instead on a characteristic time 
scale 



IV. THE USES OF THE CLOSED FORM 
SOLUTIONS OF THE FOKKER-PLANCK 
EQUATION 

The solutions presented above permit some analysis of 
several important dynamical situations. 

In many cases in field theory the harmonic solutions 
can reveal the essential physics of the (short time) evolu- 
tion. This is true for example for the motion of soft fields 
in the vicinity of a second order phase transition, where 
only the lowest relevant operators are needed, and for 
some situations when (spinodal) instabilities can develop, 
such as in the early stages of reheating after cosmological 
inflation or in the aftermath of a violent pressure quench, 
in some condensed matter experiments. 

In this section I analyze these two situations in the 
context of the solutions of section fill 



A. Critical Dynamics of the fields and the theory of 
topological defect formation at second order 
transitions 

In the vicinity of a symmetry breaking second order 
transition the physical mass squared of the fields will ap- 
proach zero, with a universal critical exponent v, 

m 2 (T)~m 2 |^r, 

where T c is the critical temperature. This can be seen 
already at first order in perturbation theory, albeit only 
with the mean- field value of v = 1/2. In 3D we have 



t eq ~ 2/77. (74) 

It is worth commenting of the form (]73). It shows that 
even in the context of a second order in time Langevin 
evolution the long wavelength modes are effectively over- 
damped in the critical domain. This is the essence of the 
perhaps more familiar TDGL evolution, which is an ex- 
pansion in the lowest number of relevant field operators 
that can only be justified rigorously in the critical domain 
of a second order transition. We see therefore that our 
Langevin description encapsulates the TDGL dynamics 
but in addition also applies more generally. 

Now, particularly in D = 3 (as compared to lower di- 
mensions) the temperature corrections to the mass are 
dominated by the small wave-length modes. Imagine 
then that the bath temperature is changed to a new value 
Tf in the vicinity of the phase transition. Over a time 
t ~ 2/77 the thermal mass will then adapt to its new small 
value m 2 (Tf). In contrast the long wave-length modes re- 
quire a much longer time to rethermalize and stay away 
from thermal equilibrium for a time t ~ rj/m 2 (Tf) >> 77. 

This imbalance is at the heart of our current under- 
standing of the dynamics of second order phase tran- 
sitions and constitutes in particular the essence of the 
Kibble-Zurek theory of topological defect formation. The 
familiar response time t(T) is here given by Eq. J7"s|). 
Incidentally the temperature dependence of ( |72|) also im- 
plies that whatever the time behavior of the external bath 
temperature the dependence of m 2 (T) on time will be the 
same, i.e. for example for a linear external drive of the 
bath temperature one also obtains m 2 (t) oc t, again on 
timescales t ~ 2/rj. 
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The consequences of this imbalance are very impor- 
tant. By crossing a second order phase transition at some 
finite rate tq imposed externally there will be a time 
in the vicinity T c when the long-wave length modes fall 
out of thermal equilibrium. Because the system always 
rethermalizes in its small scales first, soft long wavelength 
non-trivial field configurations can persist for a long time. 
In particular topological defects can be "formed" in this 
way. 

B. Spinodal instabilities 

Another interesting application of the harmonic solu- 
tions of section || is the onset of (spinodal) instabilities. 
This situation assumes an initial instability character- 
ized eg. by a coherent {i.e. quasi-spatially homogeneous) 
field, close to the origin <f> = of a double-well potential. 
In addition there may be fluctuations (thermal and/or 
quantum) about this mean field but these must be tuned 
to be small and are often taken to be Gaussian and white. 

In these circumstances m 2 is initially negative and the 
long wave-length modes k 2 < —m 2 will develop instabil- 
ities characterized by exponential growth of their ampli- 
tude at early times. In the limit of vanishing stochastic 
fields this behavior can be captured by solving the lin- 
earized Langevin equation |Ig|-|20[] , whose solutions are of 
course simple exponentials of ±iy/k 2 — m 2 t. Linear de- 
pendences of m 2 on time can also be captured in similar 
ways. 

The advantage of using the set of solutions ([47] ) in this 
context is that the noise effects are then better taken into 
account. For the linearized problem the set of eigenfunc- 
tions and eigenvalues given in pTj ) is exact. The behav- 
ior of these solutions in the overdamped case, where the 
eigenvalues are 

(k 2 — m 2 ) 

E N = N± '-, 75 

V 

is obvious: whenever k 2 — m 2 < all excited states dis- 
play exponential growth. It is worth noting that it is 
the highest lying states (i.e. those with largest N) that 
grow the fastest, leading to spectacularly out of equilib- 
rium configurations, at least relative to the initial state. 
The short wavelength modes are of course still dissipated 
away on a time scale t ~ 2/rj so these excited states are 
highly biased towards the infrared as is well known from 
direct numerical studies. 

For the second order Fokker-Planck equation the eigen- 
values are instead (pl7|). It is then clear that there will 
be growing modes as long as k 2 < — m 2 . In these cir- 
cumstances, for small but negative m 2 /r], the eigenvalues 
coincide with (|75|). For the second time derivative evo- 
lution however there will be both exponentially growing 
and decaying modes. 

The validity of the linear approximation breaks down 
when the 2-point function G(t), which controls the lead- 



ing corrections to the mass squared, becomes of order 
the negative bare mass squared. In the Fokker-Planck 
picture this is an average which correctly accounts for 
the noise effects on the short wavelength modes, which 
are quasi-Gaussian (perturbative) in any case. The value 
of G in linearized evolution can clearly be computed in 
closed form since all integrations are over polynomials of 
the fields times a Gaussian. 

Computations of eg. G(x, t) at the point where the 
linearized approximation fails have been shown to be 
useful Jig ), eg. in the computation of topological de- 
fect formation through the insertion of G(x, t) into the 
well-known Halperin formula for counting field zeros, pro- 
vided there is little energy left in the system at this point 
to cause substantial subsequent transport towards short 
wavelengths, (see Q for counterexamples). 

V. DISCUSSION AND CONCLUSIONS 

In this paper I discussed general properties of the sec- 
ond order Langevin scalar field evolutions as well as of 
some of its related models. The associated Fokker-Planck 
equations were solved in closed form for the harmonic po- 
tential and a mean-field approximation was devised for 
both second order and overdamped cases. The analo- 
gous treatment for the stochastic nonlinear Schrodinger 
equation was given in Appendix A. The latter can be 
solved numerically with an effort that is comparable to 
the mean-field Langevin equation, but has the advantage 
of already including the effects of the averages over the 
stochastic fields. 

Although limited to special cases closed form solutions 
of the Fokker-Planck equations, have been shown to be 
very useful in providing us with analytical qualitative 
(and sometimes quantitative) insights into the evolution 
of the fully non-linear field theoretical system, especially 
in determining time scales for thermalization of several 
parts of the system. Through this analysis I showed how 
fundamental ingredients of the dynamics of second order 
phase transitions are embodied in these effective models 
and how they generalize the well known TDGL evolution 
in the critical domain. 

In particular in the vicinity of T c , where m 2 (T) — ► 0, 
one obtains true time scale separation, with the small 
wavelength field modes thermalizing in a time t eq ~ if 1, 
while the long wave length modes take a time t eq ~ 
r]/m 2 (T) >> l/rj. In these circumstances, over times 
1/rj << t << rj/m 2 (T), one can take the short wave- 
length modes to be thermalized at the temperature of 
the stochastic thermal bath, and include their effects on 
the long-wavelength modes through the temperature de- 
pendence of the parameters in their potential. This treat- 
ment allows us to discuss in qualitative correct terms the 
evolution of the long wavelength modes in the critical do- 
main, the phenomenon of critical slowing down and the 
ingredients of the theory of topological defect formation 
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in these models. 
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In this appendix I repeat the analysis applied to the 
Langevin equations ([!]) to the case of a stochastic non- 
linear Schrodinger equation (sNLS), known in some in- 
stances also as the stochastic Gross-Pitaevsii equation. 
This model is important in the study of non-relativistic 
many body systems. Many examples exist that are 
thought to obey this effective dynamics, ranging from 
strong type II (or hard) superconductors and superfluid 
4 He to atomic Bosc-Einstein condensates and light prop- 
agation in non-linear media. 

The sNLS is of the form 

-(i + »?W = -VV + 5£r + (76) 

where if) = ip(x,t) is a complex held in space and time. 
The potential V is usually taken to be 

U[|Vf] = m 2 \iP(x,t)\ 2 + ±\iP(x,t)\\ (77) 

where m 2 can be positive or negative. This theory is 
essentially the non-relativistic version of the second order 
in time Langevin system for a complex scalar held. The 
theory is manifestly invariant under global phase U(l) 
transformations. Its Hamiltonian is 



h = J d D x{m\ 2 + v[\i>\ 2 }}. 



(78) 



Although the stochastic terms in Eq. ([76|) may be rep- 
resentative of additional intrinsic degrees of freedom, it 
may also be possible to actually build an experimental 
situation in which the many-body system is driven exter- 
nally and has losses to the outside so as to realize (|7^). 
In this latter picture the requirements are that the driv- 
ing held should be phase incoherent (at least over some 
short characteristic time and spatial scales) and that the 
losses would be proportional to the frequency of the ex- 
citations in the system. Perhaps the most challenging 
feature would be to find a trap or mirror whose trans- 
mission is linear in the range of frequencies expected in 
the system. 

The Fokker-Planck equation corresponding to this evo- 
lution is 



dt.P = 



rj Sip 



+ - 



SH 
> 

' 'Sip* 
1 



SI 



■i + rj Sip* 



■ 7] Sip* 

SH 

> |_ 

Sip i + 77 Sip 



(79) 



P 



where SI = 2r/T as before. The system thermalizes to its 
canonical distribution 



P =Af J Dip Dip* exp [-/3H] . 



(80) 



The set of eigenvalues and eigenvectors for the har- 
monic problem are 







E N = 2Nr, (81) 

1 + T] z 



F N = H N 



(82) 



The factor of 2 in En relative to (|5C|) accounts for 2 real 
fields in the complex quantity if), instead of one. Thus, 
the equilibration time scale is 



69 277 (fc 2 +m 2 )' 



which again displays critical slowing down for fc 2 ~ and 
m 2 (T) -> 0. 

It is clear that the overdamped limit is recovered for 
large T), as expected. We therefore see that in the over- 
damped limit all three models considered in this paper 
lead to the same characteristic time scale for equilibra- 
tion as could be expected on general grounds. The differ- 
ences arise for the short wavelength modes in the system 
and/or for small rj. 

An experimental realization of ( |76| ) would allow for 
the driving of the physical system across its phase tran- 
sition and more generally to canonical equilibrium at an 
arbitrary temperature, by changing the intensity of the 
phase incoherent driving field £. If £ had a phase coher- 
ent component, in addition to the incoherent piece, the 
drive would make the system behave like an XY mag- 
net at finite temperature in the presence of a magnetic 
field in the XY plane. This would allow for experiments 
in non-linear optics (where this effect is known as phase 
pinning) or atomic Bose-Einstein condensates to explore 
regimes similar to those realized in eg. high-Tc super- 
conductors. 
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